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The  ability  to  optically  image  with  nanometer  resolution  is  essential  for 
understanding  the  electromagnetic  properties  of  nanoscopic  objects1’3.  The  most 
common  methods  for  performing  nanoscale  imaging  are  based  on  near-field 
scanning  optical  microscopy4'8,  which  utilizes  a  tapered  probe  that  can  add 
unwanted  distortion  to  images9,10.  Such  distortion  can  be  eliminated  by  probing 
with  an  isolated,  point-like  emitter.  Here  we  present  a  method  for  performing 
nanoscale  imaging  using  a  deterministically  positioned  single  quantum  dot  (QD). 
The  QD  is  manipulated  to  nanoscale  precision  by  microfluidic  flow  control11’13, 
which  does  not  require  the  use  of  macroscopic  manipulators.  We  use  this  technique 
to  image  the  surface  plasmon  polariton  (SPP)  mode  of  a  silver  nanowire  with 
resolution  as  fine  as  10  nm  by  monitoring  the  coupling  efficiency  into  the  wire 
mode.  The  measured  field  distribution  reveals  nanoscale  features  such  as  localized 
modes  at  the  wire  ends  and  interference  in  the  propagating  SPP  waves. 

Imaging  with  a  nanoscopic  probe  presents  significant  challenges  due  to  the  lack  of 
suitable  techniques  for  nanoscale  particle  manipulation.  The  most  common  methods  for 
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particle  manipulation  are  based  on  optical  tweezers14,  but  optical  gradient  forces  scale 
with  particle  volume,  making  the  manipulation  of  nanoscopic  objects  challenging15. 
Methods  for  incorporating  nanoparticle  emitters  on  tapered-fiber  scanning  microprobe 
tips  have  been  demonstrated,  but  currently  these  techniques  have  been  limited  to  an 
imaging  resolution  of  about  100  nms.  Techniques  based  on  mechanical  placement  of 
diamond  nanocrystals  along  a  surface  have  achieved  nanoscale  positioning 
accuracy16’17,  but  so  far  have  only  provided  sparse  sampling  of  the  electromagnetic 
mode  with  too  few  data  points  to  reconstruct  an  accurate  image.  Random  diffusion  has 
also  been  used  to  decorate  plasmonic  hot  spots  with  single  molecular  emitters, 
achieving  imaging  resolution  as  fine  as  1.2  nm18.  However,  this  procedure  is  entirely 
stochastic  and  cannot  be  used  to  probe  a  desired  location  or  target  on  demand. 

In  this  work  we  demonstrate  a  method  for  imaging  with  nanometer  resolution 

using  a  single  quantum  dot  (QD)  by  utilizing  the  enhanced  electromagnetic  interactions 

between  the  QD  and  the  surface  plasmon  polariton  (SPP)  mode  of  a  silver  nanowire 

(AgNW).  QDs  are  bright  and  robust  nanoscale  emitters,  and  are  therefore  ideal 

candidates  for  probing  electromagnetic  interactions.  We  use  colloidally  suspended 

CdSe  QDs  that  are  manipulated  in  a  microfluidic  device  by  controlling  the  flow  of  the 

surrounding  liquid,  enabling  us  to  position  a  single  QD  near  the  AgNW  with  nanometric 

accuracy12.  When  a  QD  is  positioned  near  the  AgNW,  energy  from  the  photoexcited  QD 

is  non-radiatively  transferred  to  the  wire’s  SPP  mode  through  an  electric  dipole 

interaction19.  The  efficiency  of  this  process  scales  as  the  spontaneous  emission  rate  of 
the  QD  into  the  SPP  mode,  which  takes  on  the  form  TS7,  =  T0  E  (r )  ~  /  |Emax  |" ,  where 

E(r)  ~  is  the  electromagnetic  field  intensity  of  the  SPP  mode,  |Emax|~  is  the  maximum 
electric  field  intensity  of  the  mode,  and  T0  is  the  spontaneous  emission  rate  at  the  field 

maximum20.  The  electromagnetic  mode  profile  of  the  SPP  mode  can  therefore  be 
determined  by  measuring  the  rate  at  which  photons  couple  into  the  AgNW,  which  is 
monitored  by  observing  the  scattered  light  intensity  from  the  wire  end19. 
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An  optical  image  of  the  micro  fluidic  device  used  to  position  QDs  is  shown  in  Fig. 
la.  Four  tapered  microfluidic  channels  intersect  at  the  control  region  where  QD 
manipulation  occurs  (the  dashed  circle).  AgNWs  with  an  average  diameter  of  100  nm 
and  an  average  length  of  4  pm  (SEM  image  in  Fig.  lc  inset)  are  deposited  on  the 
polydimethylsiloxane  (PDMS)  surface  of  the  control  region.  The  channels  are  filled 
with  fluid  containing  QDs  that  become  confined  to  a  thin  sheath  along  the  surface  due 
to  the  fluid  chemistry13.  Within  this  sheath,  QDs  are  manipulated  to  nanometric 
precision  using  electroosmosis  with  feedback  control12’13.  The  AgNW  acts  as  an 
obstacle  for  the  QDs,  indicating  that  they  are  constrained  to  lie  within  100  nm  of  the 
surface  (the  height  of  the  AgNW).  The  device  is  mounted  on  an  inverted  microscope 
and  QDs  are  imaged  using  a  CCD  camera.  Their  positions  are  tracked  with  sub¬ 
wavelength  precision  by  fitting  the  diffraction  spot  to  a  Gaussian  point  spread  function. 
The  vision  accuracy  of  the  QD  tracking  algorithm  was  determined  to  be  9  nm  by 
tracking  an  immobilized  QD  on  a  glass  surface  (see  Supplementary  Information). 
Details  of  the  fluid  composition,  experimental  setup,  and  nanowire  synthesis  are 
provided  in  the  Methods  section. 

Fig.  lb.  illustrates  how  imaging  of  an  AgNW  is  performed.  Within  the  control 
region,  a  single  QD  is  selected  and  driven  along  a  trajectory  that  samples  the  local  field 
of  the  wire  at  a  desired  set  of  locations.  The  accuracy  within  which  QDs  can  be 
positioned  is  measured  to  be  41  nm,  ensuring  that  we  can  deterministically  probe  the 
wire  mode  at  a  desired  location  on-demand.  The  QD  trajectory  is  selected  to  maximize 
data  sampling  near  the  wire  surface  (see  Supplementary  Information  for  both 
measurement  of  the  QD  positioning  accuracy  and  selection  of  the  QD  trajectory).  Figs, 
lc-e  are  a  series  of  images  of  a  single  QD  being  moved  progressively  closer  to  an 
AgNW.  When  the  QD  is  in  close  proximity  to  the  wire,  light  scatter  is  observed  from 
the  wire  ends  (Fig.  le).  A  collection  of  such  images  can  be  used  to  map  the  SPP  mode 
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of  the  AgNW  by  monitoring  the  scattered  intensity  at  the  ends  of  the  nanowire  as  a 
function  of  the  QD  position. 

Intensity  variations  due  to  QD  blinking  as  well  as  local  field  modification  of  the 
pump  by  the  AgNW  can  complicate  interpretation  of  the  scattered  intensity.  We  account 
for  these  effects  by  normalizing  the  intensity  scattered  from  the  wire  end  by  the 
intensity  of  the  direct  QD  emission.  The  QD  emission  and  the  AgNW  scatter  are 
measured  by  summing  the  pixels  within  their  respective  windows  (Fig.  lc-e).  We  define 
the  normalized  coupled  intensity/  as: 

Z(  jNW  _  jBKG  \ 

\  n,m  ) 

T /  \  _  n,m 

1  IT  y>-  1BKG  j  (1) 

PA 

where  l^n  and  IQpDq  are  the  intensities  of  the  (n,m)th  and  (p,q)th  pixels  within  the  AgNW 

and  QD  windows  respectively,  while  IBKC  is  the  average  pixel  background  intensity. 
The  background  level  is  calculated  by  averaging  pixels  sufficiently  far  away  from  the 
AgNW  and  QD.  The  quantity  7  is  directly  proportional  to  the  electromagnetic  field 
intensity  of  the  SPP  mode. 

Fig.  2a  is  a  scatter  plot  of  7  as  a  function  of  QD  position  when  the  QD  is  near  the 
midpoint  of  the  AgNW  (video  provided  in  Supplementary  Information).  The  shaded 
region  represents  the  wire  location  as  estimated  from  light  scattered  from  the  wire  ends 
and  an  assumed  AgNW  radius  of  50  nm.  The  color  of  each  data  point  corresponds  to  7 
at  that  location,  which  is  observed  to  increase  as  the  QD  approaches  the  AgNW.  When 
the  QD  is  closest  to  the  wire,  a  maximal  7  of  0.12  is  observed.  Fig.  2b  shows  7  as  a 
function  of  radial  distance  from  the  wire  axis.  For  a  cylindrically  symmetric  AgNW,  the 

evanescent  field  from  the  surface  follows  a  Bessel  function  decay21.  A  fitting  function 
\pKn(apf  is  therefore  used,  where  K0(ap)  is  the  zeroth-order  modified  Bessel 


5 


function,  p  is  the  radial  distance  from  the  AgNW  axis,  and  a  and  P  are  fitting 
parameters  (calculated  to  be  0.01  nm"1  and  0.51,  respectively).  The  solid  blue  line  in 
Fig.  2b  is  the  Bessel-function  fit  and  the  red  dashed  line  corresponds  to  the  AgNW 
evanescent  field  as  calculated  using  finite-difference  time-domain  (FDTD)  simulation 
(see  Methods).  The  FDTD  solution  is  multiplied  by  a  normalization  constant  to  provide 
the  best  fit  to  the  data.  To  determine  the  spatial  accuracy  of  the  measured  positions,  we 
calculate  the  root  mean-squared  distance  of  the  measured  results  from  the  Bessel 
function  fit  (See  Supplement),  which  is  determined  to  be  10  nm.  This  estimate  provides 
an  upper  bound  on  the  spatial  accuracy  by  making  the  worst-case  assumption  that  all 
noise  is  due  exclusively  to  fluctuations  in  the  measured  position.  The  estimated  spatial 
accuracy  is  consistent  with  the  9  nm  vision  accuracy  calculated  previously. 

Fig.  3a  is  a  scatter  plot  of  data  recorded  when  a  QD  was  scanned  near  one  end  of 
the  AgNW  while  the  intensity  was  monitored  at  the  opposite  end.  Probing  the  AgNW 
tip  required  several  minutes  of  QD  scanning,  over  which  sample  drift  becomes 
important.  We  corrected  for  the  sample  drift  by  calculating  the  QD  position  relative  to 
the  wire  end  that  is  tracked  throughout  the  experiment  by  the  same  method  that 
monitors  the  QD  position.  The  ability  to  correct  for  drift  represents  an  important 
practical  advantage  of  this  approach,  enabling  us  to  acquire  images  over  long 
acquisition  times  without  image  distortion.  As  shown  in  Fig.  3a,  enhanced  coupling  is 
observed  near  the  end  of  the  AgNW,  which  is  attributed  to  a  localized  mode  at  the 
tip19'22.  A  maximal  I  of  0.34  is  observed,  which  is  significantly  higher  than  the 
coupling  measured  at  the  middle  of  the  wire.  Data  points  located  within  the  shaded  wire 
region  are  attributed  to  the  QD  being  pushed  slightly  on  top  of  the  wire,  which  does  not 
act  as  a  perfect  obstacle. 

The  raw  data  presented  in  Fig.  3a  can  be  used  to  construct  an  image  of  the  field 
profile  of  the  localized  mode.  The  value  of  each  pixel  in  the  image  is  found  by  taking  a 
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Gaussian-weighted  average  of  the  raw  data.  The  Gaussian  is  centered  at  the  location  of 
the  pixel  and  the  standard  deviation  is  set  to  22  nm,  corresponding  to  the  root-mean- 
square  combined  spatial  accuracy  of  the  QD  (10  nm)  and  the  tracked  AgNW  end  (20 
nm)  (see  Supplementary  Information).  The  additional  error  incurred  by  tracking  the 
wire  end  is  not  fundamental  to  the  imaging  procedure  and  could  be  largely  removed  by 
using  brighter  tracking  objects  to  monitor  the  drift.  The  resulting  two-dimensional  field 
profile  is  shown  in  Fig.  3b.  Comparison  of  the  measured  field  profile  with  the  calculated 
mode  profile  attained  from  FDTD  simulations  shows  good  agreement  between 
experiment  and  prediction  (Fig.  3c). 

In  addition  to  the  highly  localized  mode  at  the  wire  end,  Figs.  3b  and  3c  suggest 
the  presence  of  an  oscillatory  mode  structure  along  the  sides  of  the  wire.  To  examine 
this  mode  structure  in  more  detail,  we  probed  the  field  along  a  1  pm  long  region  at  one 
end  of  a  wire.  The  measurement  results  are  shown  in  Fig.  4a,  where  position  is  once 
again  plotted  relative  to  the  wire  end.  In  this  measurement  we  used  flow  to  position  a 
QD  as  close  as  possible  to  the  wire  surface.  Under  these  conditions  most  of  the  QD 
positions  were  measured  within  the  shaded  region.  These  data  points  are  again 
attributed  to  being  on  top  of  the  wire,  where  we  measured  /  values  as  high  as  0.16 
(which  is  slightly  larger  than  the  values  found  in  Fig.  2a).  This  increased  coupling 
indicates  that  the  QD  is  close  to  the  AgNW  surface. 

A  clear  periodic  pattern  is  observed  in  the  scatter  plot  as  the  QD  is  moved  along 
the  length  of  the  AgNW.  These  oscillatory  fringes  arise  from  interference  between  the 
QD  emission  propagating  in  the  forward  direction  and  reflected  from  the  wire  end1.  The 
oscillatory  pattern  is  more  readily  observed  in  the  Gaussian  image  reconstruction  in  Fig. 
4b.  Fig.  4c.  displays  the  calculated  FDTD  field  intensity  along  the  side  of  the  AgNW, 
and  exhibits  good  qualitative  agreement  with  the  measured  data.  In  Fig.  4d,  we  plot  I 
for  data  points  within  the  shaded  region  in  Fig.  4a  and  a  sinusoidal  fit.  From  the  fit  we 
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determine  wavelength  of  the  propagating  SPP  mode  to  be  368  nm.  This  value  is  slightly 
different  than  the  329  nm  calculated  with  FDTD.  This  difference  is  likely  due  to  use  of 
parameters  for  bulk  silver23  in  the  simulation.  Our  AgNWs  are  bicrystalline24  and  are 
therefore  expected  to  have  a  different  dielectric  constant  than  bulk  silver,  resulting  in  a 
different  propagation  constant. 

In  conclusion,  we  have  demonstrated  a  method  for  imaging  the  field  profile  of 
surface  plasmon  polariton  modes  with  nanometric  resolution  using  enhanced  atom-light 
interactions  with  an  isolated  single  QD.  The  demonstrated  approach  provides  a  direct 
measurement  of  the  electromagnetic  field  profile,  and  also  enables  deterministic  control 
of  these  interactions  by  precise  placement  of  the  QD  at  different  field  locations.  Fluid- 
based  methods  for  three-dimensional  tracking  and  control25  of  particles  could  also  be 
incorporated  to  perform  three-dimensional  imaging  with  nanoscale  accuracy.  In 
addition,  by  combining  this  technique  with  methods  for  selective  immobilization  of 
QDs13,  our  results  can  be  extended  from  probing  to  a  deterministic  assembly  of  active 
quantum  emitters  and  plasmonic  structures  for  development  of  new  nanoelectronic 
devices26  and  quantum  optical  circuits27. 

METHODS 

Fluid  composition.  The  control  fluid  is  composed  of  47.5%  by  volume  ethoxylated- 1 5 
trimethylolpropane  triacrylate  resin  (SR-9035,  Sartomer)28,  1.25  wt%  rheology  modifier 
(Acrysol  RM-825,  Rohm  and  Haas  Co.)29,  0.15%  of  a  zwitterionic  betaine  surfactant30, 
and  CdSe/ZnS  core-shell  QDs  (Ocean  NanoTech©,  Carboxylic  Acid,  620  nm)  in 
deionized  water.  The  triacrylate  resin  causes  the  fluid  sheath  to  form.  The  surfactant  is 
introduced  to  improve  EOF  actuation  on  the  PDMS  surface  while  the  rheology  modifier 
is  used  to  reduce  Brownian  motion  and  inhibit  unwanted  adhesion  of  the  QDs  to 


surfaces. 


Experimental  setup.  The  microfluidic  device  is  composed  of  a  molded  PDMS  cross 
pattern  placed  on  top  of  a  glass  coverslip.  The  resulting  channel  is  5  pm  in  height  with  a 
control  region  that  is  approximately  100  pm  in  diameter.  The  QDs  are  excited  by  a 
highly  defocused  532  nm  light  source  with  illumination  intensity  of  200  W/cm2,  and 
imaged  onto  an  EMCCD  camera  (Hamamatsu  C9100-13)  operating  at  20  Hz  frame  rate. 

Synthesis  of  silver  nanowires.  Silver  nanowires  were  synthesized  by  reducing  AgNCE 
with  ethylene  glycol  (EG)24.  In  a  typical  synthesis,  5  mL  of  EG  (Fisher  Scientific)  was 
placed  in  a  50  mL  round-bottom  flask  equipped  with  a  condenser,  a  thermometer  and  a 
magnetic  stirring  bar.  After  the  EG  was  heated  to  160  °C  in  an  oil  bath,  0.5  mL  of  PtCE 
(Aldrich,  1.5  xlO"4  M  in  EG)  was  injected  into  the  solution  to  synthesize  1-5  nm 
precursor  nanoparticles.  After  4  min,  2.5  mL  of  AgNCE  (Aldrich,  0.12  M  in  EG)  and  5 
mL  of  poly(vinyl  pyrrolidone),  (PVP,  Mw  ~  40,000,  City  Chemical  LLC,  0.36  M  in 
EG)  were  added  simultaneously  over  6  min  with  an  accompanying  color  change  from 
yellow  to  yellowish  gray.  When  the  AgNCE  had  been  completely  reduced  by  EG  for  60 
min,  the  solution  was  cooled  to  room  temperature.  To  remove  the  remaining  EG,  PVP 
and  silver  nanoparticles,  the  solution  was  diluted  with  ethanol  and  centrifuged  at  2000 
rpm  for  15  min.  After  repeating  several  times  until  the  supernatant  became  transparent, 
the  solution  was  redispersed  in  H2O.  The  resulting  bicrystalline24  silver  nanowires  are 
100  to  200  nm  in  diameter  and  3  to  10  pm  in  length,  as  determined  by  scanning  electron 
microscopy. 

Finite  different  time  domain  simulation.  Simulations  were  performed  using  the 
Lumerical  FDTD  Solutions  software  package  (http://www.lumerical.com).  A  4  pm  long 
silver  nanowire  with  100  nm  diameter  was  simulated  surrounded  by  a  background  index 
of  1 .4,  which  corresponds  to  the  index  of  refraction  for  PDMS.  The  QD  fluid  was 
assumed  to  have  the  same  index.  The  actual  index  of  the  fluid  lies  somewhere  between 
water  (1.3)  and  the  index  of  the  resin  (1 .47). 
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Figure  1  |  QD  near-field  probing,  a.  Optical  image  of  the  microfluidic  crossed- 
channel  device.  Flow  in  the  center  control  region  (dashed  circle)  is  manipulated 
in  two  dimensions  by  4  external  electrodes  (not  shown),  b.  Illustration  of 
positioning  and  imaging  technique.  A  single  QD  is  driven  along  a  trajectory 
close  to  the  wire  by  flow  control  and  the  coupling  between  the  QD  and  AgNW  is 
measured  by  the  scattered  intensity  from  the  wire  ends.  The  inset  shows  an 
SEM  image  of  a  typical  AgNW  used  in  our  experiments,  c-e.  A  series  of 
emission  images  showing  coupling  of  the  QD  to  the  AgNW  as  the  QD  is  moved 
closer  to  the  wire.  Intensities  are  plotted  on  a  log  scale.  Red  (blue)  boxes  signify 
the  image  integration  region  used  to  calculate  the  intensities  and  positions  of 
the  QD  (wire  scattering  points).  The  measured  location  of  the  QD  is  labeled  with 
a  red  star  and  the  axis  of  the  AgNW  is  labeled  with  a  blue  dashed  line. 
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Figure  2  |  Probing  the  evanescent  field  of  an  AgNW.  a.  I  as  a  function  of 
position  near  the  middle  of  the  wire.  The  shaded  area  indicates  the  physical 
extent  of  the  AgNW.  b.  I  as  a  function  of  distance  from  wire  axis  using  data 
from  panel  a.  The  blue  line  indicates  the  best  fit  to  a  modified  Bessel  function 
The  red  line  is  an  FDTD  simulation  of  the  AgNW  evanescent  field.  The 
simulation  result  was  fit  to  the  data  using  an  overall  scaling  factor. 
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Figure  3  |  Evanescent  field  at  the  tip  of  an  AgNW.  a.  Scatter  plot  of 
measured  positions  and  intensities  near  the  end  of  the  AgNW.  The  color  of 
each  data  point  corresponds  to  7 .  The  shaded  area  indicates  the  position  of  the 
wire.  b.  Reconstructed  image  using  a  Gaussian  weighted  average.  The  image 
intensity  is  normalized  by  its  maximum,  c.  FDTD  simulation  of  the  AgNW  mode 
profile  showing  an  enhancement  at  the  tip  (also  normalized  by  its  maximum). 
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Figure  4  |  Probing  the  mode  structure  of  an  AgNW.  a.  Scatter  plot  of  the 
measured  QD  positions  near  the  end  of  the  wire.  The  color  of  each  data  point 
corresponds  to  /  ,  measured  at  each  location.  The  shaded  area  indicates  the 
position  of  the  AgNW.  b.  Reconstructed  image  using  a  Gaussian  weighted 
average.  Image  intensity  is  normalized  by  its  maximum,  c.  FDTD  simulation  of 
the  field  intensity  standing-wave  pattern  along  the  side  of  the  AgNW  (also 
normalized  by  its  maximum),  d.  Circles  indicate  the  measured  value  7  as  a 
function  of  position  along  the  wire.  The  solid  line  represents  a  sinusoidal  fit  to 
the  data. 
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